Atomic Fermi gas in the trimerised kagome lattice at the fiUing 2/3 
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We study low temperature properties of a spinless interacting Fermi gas in the trimerised kagome 
lattice. The case of two fermions per trimer is described by a quantum spin 1/2 model on the 
triangular lattice with couplings depending on the bond directions. Using exact diagonalisations we 
show that the system exhibits non-standard properties of a quantum spin-liquid crystal, combining 
a planar antiferromagnetic order with an exceptionally large number of low energy excitations. 

PACS numbers: 03.75.Ss,05.30.Fk 
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One of the most fascinating recent trends in the physics 
of ultracold gases concerns atomic gases in optical lat- 
tices, where strongly correlated systems may be realized. 
Such systems offer an "atomic Hubbard toolbox" 1] to 
simulate various sorts of Hubbard models, and to study 
phenomena known in condensed matter physics in an un- 
precedently controlled manner. To name just few exam- 
ples, atomic lattice gases may serve to study various spin 
models 0, to simulate high Tc superconductivity to 
investigate a variety of quantum disordered systems 0, 
or to process quantum information Q- Seminal exper- 
iments of Ref. have stimulated a great interest in 
experimental studies of atomic lattice gases (cf. 

Particularly fascinating in this context is the possi- 
bility of studying quantum frustrated antiferromagnets, 
which lie at the heart of modern quantum magnetism 
0. Recently we have proposed how to create ideal 
and trimerised kagome optical lattices, and have stud- 
ied physics of various quantum gases in such lattices . 
A Fermi-Fermi mixture with half filling for both species 
in the limit of strong interspecies coupling behaves in 
the ideal kagome lattice as a spin 1/2 Heisenberg an- 
tiferromagnet (KAF). Such a system (having so far no 
experimental realization among solid state systems) is 
a paradigmatic, although not yet fully understood (cf. 
01) example of a quantum spin liquid of type H 0. 

In Ref. we have also discussed briefly the case of 
interacting spinless Fermi gas in the trimerised kagome 
lattice at filling 2/3 (2 atoms per trimer). Such a system 
behaves as a quantum magnet on the triangular lattice 
with couplings that depend on bond directions, and is 
particularly interesting since: i) it describes the physics 
of the trimerised KAF in the plateau region at one third 
of the saturation magnetisation ; ii) it has itself fas- 
cinating properties, expected to be generic for "multi"- 
merised frustrated systems; iii) it is a paradigmatic Fermi 
system to study in trimerised lattices; iv) it is experimen- 
tally feasible. In this Letter we study low temperature 
T physics of this system using exact diagonalisations of 
the Hamiltonian for 12 • • • 24 spins. We show that for ef- 



fectively ferromagnetic couplings the system exhibits the 
non-standard properties of a quantum spin-liquid crystal, 
combining planar antiferromagnetic order with an excep- 
tionally large number of low energy excitations, and a 
small (if any) gap. For the effectively antiferromagnetic 
coupling the quantum results agree very well with the 
classical results indicating antiferromagnetic planar or- 
der and a gapped spectrum. 

The experimental realization of the considered system 
requires a creation of trimerised kagome lattice, using su- 
perlattice techniques as shown in Ref. P|. The spinless 
interacting Fermi gas can then be formed, for instance, in 
a Bose-Fcrmi mixture, in the strong coupling limit, when 
bosons form a Mott insulator (MI), and fermions together 
with 0, 1, ... bosons (bosonic holes) form fermionic com- 
posites Jl^l- Alternatively, one could use a gas of po- 
larised ultracold dipolar fermions, that interact via re- 
pulsive dipolar potential. 

The spinless interacting Fermi gas in the trimerised 
kagome lattice is described by the extended Fermi- 
Hubbard Hamiltonian i/pH = - Y.{ab){^abflh + h.c.) -|- 
X](ab) UatnaTib, whcrc a = {a, i} with a referring to intra- 
trimer indices and i numbering the trimers. The tab 
and Uab take the values t and U for intra-, and t' and 
U' for inter-trimer hopping, Ua — flfa, and fa is the 
fermionic annihilation operator. The sites in each trimer 
are enumerated as in Fig. We denote the 3 differ- 
ent intra-trimer modes by = (/i,,, -I- f2.i + f3.i)/V3 
(zero momentum mode), and fjf' = (/i,i -f ^±/2,i + 
z±f3,i)/ (left and right chirality modes), where z± — 
exp(±27ri/3). 

In the limit of weak coupling between the trimers of 
the original kagome lattice, the problem of two fermions 
per trimer (filling 2/3) becomes equivalent to a quantum 
magnet on a triangular lattice with couplings that depend 
on the bond directions as described by the Hamiltonian 
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where TV denotes number of trimers, J = 4U' /9, and 
the nearest neighbours are enumerated as in Fig. In 
Eq. Q we have Si{(t)) = cos((/))si'^ + sm{(/))sy\ where 
the spin-1/2 operators Sx\ Sy"^ are defined as: si*^ = 

(/fVi^^ + /I'^Vf )/2, 4'^ = -^(/f Vl^^ - )/2. 

The angles (/> are: = </)i^6 = 0, (l>i^2 = (f>i^3 = 

27r/3, 0i_-»4 = (?!)i-.5_= -27r/3, = (/ii_^2 = -27r/3, 

(/>i— >3 = 0i— >4 = 0, (f>i^5 — 4>i^e — 27r/3. This Hamilto- 
nian has previously appeared in the context of a block- 
spin approach to the Heisenberg K AF [l3[ . The main 
purpose of that approach has been to find the origin of 
the exponentially large number of low-lying singlets that 
had been found in numerical studies of the kagome an- 
tiferromagnet HI [3. From Refs. [H [13 it also fol- 
lows that Htrimer describcs the physics of the trimcrised 
KAF in a magnetic field that drives this system into the 
plateau region at 1/3 of the saturation magnetisation. 
We stress that the Hamiltonian Htrimer to be studied in 
this Letter describes a physically feasible situation. 

Let us begin by discussing the classical theory of the 
model which describes the large spin limit. In ad- 
dition to being translationally invariant, the model of 
Eq. Q is invariant under the point group of order 6, 
Zq = Z3 ■ Z2, where the generator of Z3 (order 3) is the 
combined rotation of the lattice by the angle 47r/3, and of 
the spins by the angle 27r/3, while the generator of Z2 (or- 
der 2) is the spin inversion in the x — y plane. We remark 
that our model possesses no continuous spin rotational 
symmetry. There exist three ordered classical states with 
small unit cells that are compatible with this point-group 
symmetry of the model: a ferromagnetic state, and two 
120° Neel type structures with left (Fig. QJd) and right 
(Fig. 12^^) chiralities. The energies per site of these states 
are i^^f™ = = -35V/4 and i^^f^L = 35V/2, 

where the subscripts "right" and "left" refer to chirali- 
ties. Hence, for J < the state with left-handed chirality 
will be the ground state (GS). For J > the situation is 
more complex: the states with right-handed chirality and 
the ferromagnetic state are degenerate ground states. 

To understand further the nature of classical ground 
states we have done a numerical analysis of the 12- 
spin cell by fixing the direction of every spin to mr/S 
(n = - • -5), and checking the energies of the resulting 
6^^ configurations. This analysis has revealed that for 
J < there are 6 ground states {Zq symmetry of (^) 
each of them exhibiting the left chirality Neel order. For 
J > the results are dramatically different: there are 
240 degenerate classical GSs in this case among them 6 
pure right chirality Neel states and 6 purely ferromag- 
netic states. For an illustration, we show in Fig. |3| two 
ordered GSs with very large unit cells (Figs. Ot), EJi) to- 
gether with their parent states (Figs. EjiiEt). As will be 
seen below the large number of degenerate classical GSs 
finds its analogue in a large density of low-lying excita- 
tions of the quantum version of Eq. Q • 
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FIG. 2: (colour online), (a) classical 120° state with right 
chirality. Dots show s^ ^ s^'^^ -I- Sy 'sy^"' , where |s | = 1/2 p^: 
(b) spin-spin correlations, {sx^ s^J''^^ + s^'sj,^"''). The upper 
[lower] set of values corresponds to fcT = [kT = 10~^J/2]. 
In both plots A'^ = 21 and i = 10 at the central site. 



To describe the physics of spinless fermions on a 
trimerised optical kagome lattice at filling 2/3 we need to 
consider the model, Eq. for spin 1/2, i.e. in the ex- 
treme quantum limit. Questions to be answered for this 
case are: (z) Is the GS of the model Eq. ^ an ordered 
state, or is it a spin liquid either of type I, i.e. a state 
without broken symmetry, with exponentially fast decay- 
ing spin-pair correlations and a gap to the first excitation, 
or of type II, i.e. a kagome-like GS again without bro- 
ken symmetry, with extremely short ranged correlations, 
but with a dense spectrum of excitations adjacent to the 
GS; (ii) What are the thermal properties of our system? 
After all, the model can only be realized at finite, albeit 
low temperatures. 

To answer the above questions we have performed ex- 
act diagonalisation of the Hamiltonian for N= 12, 15, 
18, 21, and 24 spins using the ARPACK routines [lij . 
To simplify the calculations we block-diagonalised the 
Hamiltonian by exploiting its translational symme- 
tries thereby reducing the dimensions of the matrices that 
had to be diagonahsed from 2^ x 2^ to w 2^/7V x 2^/iV. 
Despite all these efforts, studies of larger systems require 
the use of massive computer resources. Fortunately, the 
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FIG. 3: (colour online) (a) Right chirality Neel configuration; 
(b) localised defect in configuration (a) ; (c) ferromagnetic con- 
figuration; (d) line defect in configuration (c). Open arrows 
present spins determined by the boundary conditions for the 
12-spin cell. Defects are marked by red dashed contours. 
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TABLE L Spin correlations for J > 0. For every distance 
from a reference centre (1, . . . , 3 in lattice units) the left most 
number is a planar spin-spin correlation, while the right most 
one is the absolute value of that correlation divided by the 
nearest neighbour spin correlation in the considered system. 



results for 21 and 24 spins show qualitative and quantita- 
tive resemblance, and we regard them as representative 
for larger systems. 

Our findings are presented in Figs. \^ and 01 and Ta- 
bles I and II. For J > 0, in contrast to the classical result, 
the ground state exhibits the 120° Neel order with right 
chirality This is illustrated in Fig. I^b, where the 

planar spin-spin correlations are presented. Direct com- 
parison with the correlations of the classical state. Fig. 
Et, shows that the exact quantum correlations, although 
smaller, have the same order of magnitude and sign as 
the classical ones. Especially, the relative values of cor- 
relations compare nicely to the classical result: Table I 
summarises the results for iV = 21 and TV = 24. Amaz- 
ingly, the 120° Neel order survives at finite temperatures, 
as is indicated by the results obtained for kT = 10~^ J/2, 
see Fig. Eb. At such temperatures about 800 low energy 
eigenstates contribute to the correlations. For smaller 
systems, < 21 (J > 0), finite size effects affect the 
spin correlations strongly. Nevertheless, the ground state 
energy per spin can be reliably extracted from the data 
for N > 12, resulting in -0.2175 J - 0.0755 J/N. 

The selection of a GS with 120° planar Neel order from 
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FIG. 4: (colour online). Number of states in an energy in- 
terval above the GS; black (solid), red (dotted) and green 
(dashed) lines for N = 24, 21 and 18, respectively. 



the large manifold of classical GSs by quantum effects im- 
plies the breaking of the translational and of the point 
group of our model Eq. but there is no continuous 
symmetry that the ordered GS could break. Therefore, 
the standard expectation would be that the excitations 
have a gap of order of J. Instead, we find that the system 
has an exceptionally large number of low energy excita- 
tions (see Fig.0J. For instance, for A^ = 21 in the energy 
interval O.lJ/2 there are about 800 excited states. Most 
of them support the spin order of the GS so that this 
order persists at finite temperatures. 

The analysis of the results for different A^'s are com- 
patible with an exponential increase of the number of 
low-energy states with the system size A'^ similarly as in 
the case of the S" = 1/2 KAF [H [3. For the KAF 
Mila has been able to explain this high density of low en- 
ergy states by associating them approximately to dimer 
coverings of an effective triangular lattice with uncorre- 
lated products of nearest-neighbour pair states 0| . His 
method fails here, because the low-lying states of our 
model must certainly be highly correlated. On account 
of the breaking of the discrete symmetries of our model 
by the Neel order, one expects the ground state of the in- 
finite system to be sixfold degenerate. For finite systems 
this degeneracy is lifted. Nevertheless, we expect to find 
six low-lying states in the gap below the lowest excited 
state. In view of this scenario the inspection of the lower 
panel of Fig. ^ suggests that the gap, if any, is smaller 
than 10~^ J/2. The appearance of this very small energy 
scale is completely unexpected and puzzling. Obviously, 
the answer to the questions (?) and (ii) above, is that the 
GS is ordered, and that the order survives at low T. The 
smallness of the gap and the large density of low-energy 
states, however, resemble very much the behaviour of 
a quantum spin liquid of type II. For these reasons we 
propose to term our system a quantum spin-liquid crys- 
tal. We remark in this context, that the specific heat of 
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TABLE II; Spin-spin planar correlations for J < as a func- 
tion of distance 1, . . . , V7 in lattice units. 



iV = 21 system exhibits a peak at fcT 7 • 10"^ J/2. 

The above results for J > contrast dramatically with 
those for J < 0, summarised in Table II. In the latter 
case we deal with the standard quantum antiferromag- 
net with 120° Neel order and left chirality I?] (Fig. [Hd). 
The spectrum is gapped, and the classical spin-spin cor- 
relations approximate the quantum correlations remark- 
ably well. The gap is of the order |J|/2 in this case 
(A^ = 12,18,21), meaning that there are at most a few 
states with energies substantially below \ J\/2 for J < 0, 
as opposed to the huge number for J > (Fig.^. 

The observation of physics described in this Letter re- 
quires achieving low, but not unrealistic T ~ lOnK to 
lOOnK (c.f. Ref. ^3)- In experiments could vary from 
~ 20 to ~ 1000. The low energy states may be prepared 
by employing adiabatic changes of the degree of trimeri- 
sation of the lattice. For instance one can start with a 
completely trimerised lattice; the filling v = 2/3 may be 
achieved then by starting with v = I, and eliminating 
1 atom per trimer using, for instance, laser excitations. 
One can then increase t and U slowly, on the time scale 
slower than the final 1/ J (~seconds) . Alternatively, one 
could start with i/ ~ 2/3 in the moderately trimerised 
regime. As in Ref. 0, the inhomogeneity of the lattice 
due to the trapping potential, would then allow to achieve 
the Mott state with i/ = 2/3 per trimer in the centre 
of the trap. Nearly perfect 2/3 filling can be achieved 
by loading a BEC of molecules formed by 2 fermions 
into a triangular lattice, generating an MI state, adia- 
batically transforming the lattice to a trimerised kagome 
one, "dissociating" the molecules by changing the scat- 
tering length to negative values, and by finally optically 
pumping the atoms into a single internal state. Prepar- 
ing = 2/3 might involve undesired heating (due to op- 
tical pumping), which can be overcome by using laser, 
or phonon cooling afterwards (cf. |13|)- Note that the 
imperfections of v can be described by a — J" -kind of 
model, and are of interest themselves. 

After state preparation it should then be possible to 
measure the energy of the system simply by opening the 
lattice; by repeated measurement of the energy E{T) at 
(definite) finite temperatures one would get in this way 
an access to the density of modes, i.e. could compare 
the results with Fig. ^ From such measurements one 
could infer about the existence of the gap i^gap, since if 
i?gap is large enough, E(T) becomes T-independent for 
kT < i?gap. Various other correlations could be mea- 
sured using the methods proposed in Ref. . In order 



to measure planar spin correlations, one has, however, 
to lift the degeneracy of the f± modes, e.g. by slightly 
modifying the intensity of one of the superlattices form- 
ing the trimerised lattice. This should be done on a 
time scale faster than the characteristic time scales of 
other interactions, so that the state of the system would 
not change during the measurement. In such a case one 
can use far off resonant Raman scattering (or scattering 
of matter waves) to measure the dynamic structure fac- 
tor, which is proportional to the spatio-temporal Fourier 
transform of the density-density correlations. At frequen- 
cies close to the two photon Raman resonance between 
the f± modes, only /_|_-/_ transitions contribute to the 
signal, and hence such measurement yields the desired in- 
formation about the correlations (/^''V-''/- ^^^d 
the spin correlations of Fig. 2. 
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